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STABILITY ESTIMATES IN INVERSE SCATTERING 

A. G. Ramm, Mathematics Department, 
Kansas State University, Manhattan, KS 66506-2602, USA 
Abstract: An algorithm is given for calculating the solution to the 3D inverse scattering problem with 
noisy discrete fixed energy data. The error estimates for the calculated solution are derived. The methods 
developed are of general nature and can be used in many applications: in nondestructive evaluation and 
remote sensing, in geophysical exploration, medical diagnostics and technology. 

I. Introduction 

Let q E Q {q: q{x) = for \x\ > a, x G R^, q{x) = q{x), q G L°°}, a > is a constant and 
the bar stands for complex conjugate. Consider the equation 

IqU ■.= [W^ + l-q{x)]u{x,a)^0 in R^, a e (1.1) 

where S'^ is the unit sphere in and 

u{x, a) — exp{ia ■ x) + exp{ir)A{a' , a) + o{r^^), r = |a;| ^ oo, a' — x/r (1-2) 

The existence and uniqueness of the solution to (1), (2) are well known and are easy to prove (see e.g. [1,2]). 
The coefficient A{a' ,a) is called the scattering amplitude at the fixed energy (/c^ = 1). Its properties are 
discussed e.g. in [2], pp. 219-246. The inverse scattering problem (ISP) consists of finding q (z Q from the 
knowledge of A{a',a) for all a', a G S^. Uniqueness of the solution to ISP is proved in [3,4], (see also [2]) 
even in the case when A{a' , a) is known for a G 5^, a' G S"!, where S^, j = 1,2, are arbitrary small open 
sets in S*^. 

The uniqueness result holds with the same proof as in [4] for the problem 

[V^ + 1 - qo{x) - q{x)]u = in i?3 and (1.2) holds (1.1') 

where qo^x) is a known short-range potential, e.g. |9o(a;)| < c(l -|- |a::|)~^^^, e > and q{x) G Q (this was 
noticed by R. Weder, in a preprint, 1991). The aim of this paper is to complete the study of the stability of 
the solution to ISP with noisy data which the author carried over in the papers [4-18] and in the monograph 
[1]. The statement of this problem (ISPs) is as follows. 

Let a function As {a', a) be given which is not assumed to be a scattering amplitude corresponding to 
a potential and let 

sup \As{a',a) - A{a',a)\ < 6 (1.3) 

The inverse scattering problem with noisy data is: 

(ISP5): Given As{a' , a) for all a', a G S"^ find qs{x) such that 

sup \qs - q{X)\ < Tj{S) ^ as S (1.4) 

where 

q{X) := / exp(— zA • x)q{x) dx. 

In other words, we want to calculate a stable approximation to q{X) and estimate the rate of decay of ri{S) 
as (5 — > 0. To get such an estimate uniform in q belonging to a compact subset of Q we assume that 

qeBc:= {q: Ml^ + ||Vg||L- < C} (1.5) 



Key words and phrases: Inverse scattering, stability, noisy data 
AMS Mathematics subject classification 35R30 (1991 Revision) 



1 



The ISP,5 is of considerable practical interest: it is a theoretical basis for many problems in nondestructive 
evaluation and remote sensing, including, to mention a few, problems of geophysical exploration, medical 
diagnostics, technology, etc. To the author's knowledge, there were no results concerning three-dimensional 
ISP^ except those given in the works [1,4-18] 

The basic results of this paper are new. They are formulated as Theorems A, B and C in Section 
II. Theorem A gives a rigorous inversion procedure for ISP with exact data and an error estimate for this 
procedure. Theorem B gives a rigorous contruction of the solution to ISPa, and an estimate for the difference 
between this solution and the unknown potential. Theorem B is of practical interest and can be used in 
designing a numerical code for solving ISP5. An alternative method is given in Section VII. 3. Theorem C 
gives an estimate for the difference of the Fourier transforms of two potentials if an estimate for the difference 
of the corresponding scattering amplitudes is given at a fixed energy. 

There are other new results in this paper. They are formulated as lemmas used in the proof of the 
basic results. Although the ideas and the techniques from works [4-18] are used in this paper there are 
many improvements, simplifications in the arguments and the paper is essentially self-contained: no prior 
knowledge of the author's work is assumed. The results may be of interest to mathematicians, engineers, 
physicists and numerical analysts interested in ISP5. 

The paper is organized as follows: Section II contains the statement of the basic results and a number 
lemmas. Section III contains auxiliary results and proofs of Theorem A and lemmas 10 through 17. Section 
IV contains proofs of lemmas 1 through 9 except lemma 6, which is proved in Section V. In Section VI proofs 
of Theorems B and C are given. In Section VII some numerical aspects of the ISP^ are discussed and a 
summary of the numerical approach is given. In addition, an alternative numerical method for solving ISPg 
is discussed and an error estimate for this method is proved. 

II Formulation of the results 

1. First, assume that the exact data are given: 

00 

A{a',a) =J2Ma)yi{'^'), Ma) ■= / A{a' ,a)Yi{a') da' (2.1) 

Here and below the summation in I denotes Yl%^oY^m=-V -^(■{(^) = -Aimio), ^^(q;) = ^^^(a), 
< m < ^, 



1 



(2£+l)(^-|m|: 



1/2 



-P^,|m|(cos^?) exp(im^), (2.2) 



(^+|m|)! 

where {"&, (j)) are the angles in the spherical coordinates determining the unit vector a e 5^: 

Q!i = sin 1? cos 0, a2 = sin 1? sin 0, as = cosi? (2.3) 

Define 

M:={6':6'eC^, 61- 61 = 1}, (2.4) 



6'-w = 6'i -wi + eis-wa + eg-wa, |6i| = (0 • e)^/^ (2.5) 

liO < d < TT, < (j) < 27r, then the vector a = (ai,a2,Q;3) runs through S^. If run through the 
complex plane C then the corresponding a runs through M' c M. The subset M' contains all the vectors 
of M except the ones of the form (vi,V2., 1) and {vi,V2, —1), where v := (wi,f2) € are vectors with the 
property v ■ v = 0, and v ^ 0. From (2.2) it follows that Yem{c() is defined for all a G M': it is sufficient to 
find for a G M' the corresponding complex numbers and (p such that formulas (2.3) hold and calculate the 
right-hand side of (2.2) for these "& and ^. Obviously ex.p{im(f)) is defined for complex (j), and the function 

P,„,(c».)Hsi„.)Mf!^M (-) 

is defined for complex 1?. Let 

Ba := {a;: x e R^, \x\ < a}, Qa ■■= {q--q = q, q = for \x\ > a, q G L^{Ba)}. 
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sup |AU«)l<ca(^j ^ (2.7) 



aeS 

-e<m<e 



|r,(0)| < ^H^M, ,>o, n:=\lm0\ (2.8) 



Lemma 1 : If q €: Qa then 



Lemma 2: If 6 € M' then 



where r > is an arbitrary fixed number and ji{r) is the spherical Bessel function, 
je{r) ■= ['^/i'^r)V^^Je+i{r), J^+i(r") is the usual Bessel function. 
It is well known that 

2^ + 1 

je{r) = ^(-^] ' , ^ [1 + 0(1)] as Z ^ oo (2.9) 

uniformly in r e [0, a] for any fixed a > 0. Prom (2.7)-(2.9) it follows that the series 

00 

A{e',a) = ^Ae{a)Ye{e'), O'&M, a e (2.10) 
e=o 

converges absolutely and uniformly on M x 5^, where M C M' is an arbitrary compact subset in M'. Indeed, 
choose r > a in (2.8). Then the series (2.10) is majorized by the convergent series 

00 2i + l 

c exp(Kr) ^ y2ZTl (-) ' (2.11) 

Note that given the data A{a',a) \/a' ,a G S"^ one finds the Fourier coefficients Ai{a) by formula (2.1) and 
defines A(6»',a) W e M, a G 5^ by formula (2.10). 

2. We now pass to the description of the inversion formula for exact data. Define 

p{u) ■.= exp{-ie-x) j u{x,a)i^{a)da-l, 9 G M' , u G L^{S^) (2.12) 
Consider the variational problem 

||/>(i/)||:=inf:=d(0) (2.13) 
where p{p) is a function of x G and 9 G M', 

IIpII := l|p||L2(B,\Bai)' a<ai<b (2.14) 

and ai,b arc arbitrary numbers subject to the restriction (2.14). Note that in the annulus ai < |a;| < 6 the 
scattering solution u{x, a) is defined explicitly through the data. Namely the following simple lemma holds: 
Lemma 3: One has 

00 

u{x, a) = exp(iQ • a;) + Ai{a)Y^{a')h(,{r), r > a, r = \x\, a' = x/r (2-15) 

£=0 

where h({r) are the spherical Hankel functions normalized so that hi{r) ~ exp(ir) as r — > 00. 
Lemma 4: If q G Qa then 

d{9) < c\9\-'^ , |6'|>1, 9gM' (2.16) 
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Here and below c stand for various positive constants which do not depend on 6 and 5. These constants 
depend on the norm := ||5||l2(b„) for q € Q, and for q € Be they depend on C where C is defined in 
(1.5). The notation \6\ 1 means that 6 is sufficiently large. 

Denote by v'{a,6) an arbitrary function in L'^{S^) which satisfies the inequality 

\\p{u{a,e))\\<{c+m-\ 1^1 ^00 (2.17) 

It follows from Lemma 4 that such a function can be calculated by solving problem (2.13). Define for this 
iy{a, 9) the quantity 

q:= -An j A{9',a)i^{a,9)da, 9' ,6 & M' (2.18) 

Fix an arbitrary large number Aq > 0, take any A e i?^, |A| < Aq, and pick any 6' , 6 such that 

e'-9 = \ 9',9eM, \9\^oo (2.19) 
Lemma 5: For any A G i?", n > 3 there exist 9,9' G M satisfying (2.19). 

Remark: For n = 2 this is not true as one can easily check. If A3 ^ ±1 then one can find 9, 9' G M' 
satisfying (2.19). In what follows we may assume without loss of generality that A3 ^ ±1 and use M' 

everywhere. Indeed, we can find the Fourier transform of the potential, q{\) for any A with A3 7^ ±1 and 
since q is continuous, in fact analytic, with respect to A, the value of g at A = (0,0, ±1) is determined 
uniquely by the continuity of q. 
Let 

g(A) := j q{x) ex.p{-iX ■ x) dx (2.20) 

Ba 

3. We are ready to formulate the first theorem. 
Theorem A: If q e Q, q is defined by (2.18) and (2.17), (2.19) liold, then 

sup \qiX) - q\ < c\9\-\ \9\ ^ ^ (2.21) 

|A|<Ao 

The constant c in (2.21 ) can be chosen uniformly in q G Be where Be is defined in (1.5). 
Corollctry: The following inversion formula holds 

q{X) = -4tt lim / A{e' , a)iy{a, 6) da, 6' - 6 = X, 6' ,6 & M' (2.22) 

Estimate (2.21) is the error estimate for the inversion formula (2.18) for the exact data A{a',a). Note 
that as the data we actually use the coeflacients Ai{a), £ > 0. In the proof of Theorem A the following 
major lemma is used: 

Lemma 6: If p is defined by (2.12) and \9\ > 1, 6» G M, then 

Ma < c{\\p\\ + \9r), II • lU := II • h^B^) (2.23) 



4. We will turn now to the description of the inversion formulas for noisy data As(a',a) given for all 
a',a G S"^. It will become clear later that we actually can use only a discrete subset of noisy data As{aj,ap) 
for some a'j,ap G 5^. Let [x] denote the integer nearest to the real number x. Define 
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As{9',a) := ^ Aseia)Ye{e'), Ase{a) := / As{a' ,a)Ye{a') da' (2.25) 



N(S) 



us{x, a) := exp{ia ■ x) + ^ Ase{a)Ye{a')he{r) (2.26) 

P5{v) := exp{-i0 ■ x) j us{x, a)u{a) da-1, 6 e M' (2.27) 

At(5) = exp[-jN{d)], 7 := In ^ > 0, (2.28) 

a{u) := Ml2^s-) (2.29) 

Lemma 7: One has, with k := \Im 6\, 

\\ps{i^)\\ < \\p{i^)\\ + ca{u) exp{Kb)n{6), (2.30) 

\\p{i^)\\ < \\ps{i^)\\ + ca{u) exp(K&)/x((5). (2.31) 

Consider the variational problem 

|6l| = sup := e(<^), \\ps{v)\\+a{v)exp{Kb)^l{5) <c\e\-^, € M' (2.32) 

where c > is a sufficiently large fixed constant, the supremum is taken over u e LF'{S'^) and 6 e M' , 5 > Q 
being fixed. 
Lemma 8: 

Q{5) ^ 00 as (5 ^ (2.33) 

See also formula (4.39) for a stronger result. 

Let vs{ct) G L'^{S'^) be any function such that 

\\p6{y)\\+a{v,)exv{Kh),xr{5) < c\6{5)\-\ \e{5)\ > ^ (2.34) 

where k = n{5) := |Im 9{5)\. One can calculate such vsict) and 9{6) by solving problem (2.32). Find 
e'{5) e M' such that (2.19) holds with an arbitrary fixed A e B?, |A| < Aq. Define 



:=-47r j As{e'{5),a)y5{a)da (2.35) 



Our main result is an error estimate for the inversion formula (2.35) for noisy data. We assume that (1.3) 
holds and q ^ Q where Q is defined in Section I. 
Theorem B: If q e Q and (1.3) holds, then 

sup \q{X) - qs\ < c\e{S)\-' (2.36) 

|A|<Ao 

where qs and q{X) are deGned in (2.35) and (2.20), 
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the constant c in (2.36) docs not depend on 5 and can he chosen uniformly for q e Be- 

5. If one applies a quadrature formula to (2.35) one obtains 

n 

qs « -ATrJ2M0'{S), aj)iys{aj)cj (2.38) 

where aj arc the nodes and cj are the coefficients of a quadrature formula. Formula (2.38) uses discrete 
noisy data. Using the definition (2.25) and applying a quadrature formula again one has 

^e{aj) ~ Ma'p, aj)W^c'p (2.39) 
p=i 

Combining (2.35), (2.25), (2.38) and (2.39) one gets 

n N(S) m 

q « -A^Y.^s{aj)cj ^ Y,{0\5))Y,As{a'^,ai)Y^% (2.39') 

j=l t=0 p=l 

This formula uses the discrete noisy data As{a'p,aj), I < p < ni, 1 < j < n. In principle one can estimate 
the error of the quadrature formulas (2.37)-(2.39') in terms of some bounds on the derivatives of As{6'{6), a) 
and vsict), but we do not go into detail. 

6. Let us consider the following question. Suppose qj & Q, j = 1, 2, Aj{a' , a) is the scattering amplitude 
corresponding to the potential qj. The question is: how does one estimate p := qi — q^ in terms of 
A:= Ai — A2. The following lemma is useful. 

Lemma 9: One has 

— 47ryl(a', a) = p{x)ui{x,a)u2ix,—a')dx (2.40) 



where Uj{x, a) is the scattering solution corresponding to the potential qj{x), j = 1, 2. 

The answer to the above question is given in the following theorem. 
Theorem C: If 

sup \Ai{a' ,a)- A2{a' ,a)\<5 (2.41) 

then 

sup \m\ < (2.42) 
|A|<Ao |ln<)| 

where c=const> can be chosen uniformly for q G Be- 
lli. Auxiliary Results and Proofs of Theorem A and Lemmas 10—17 

1 . In this section we give some auxiliary results which are of independent interest and also help to prove 
the results formulated in Theorems A-C. 

Let us start with an estimate of the quantity q{X). The starting point is the well-known formula 

-4:7rA{9',a) = j exp{-i0' ■ x)u{x,a)q{x) dx (3.1) 

Multiply (3.1) by i'{a) G L^{S^) and integrate over S^. Write the result as 

-47r j A{d', a)u{a) da = j exp[-i(6l' -e)-x][l + p{u)]q{x) dx (3.2) 

Here p{u) is defined in (2.12), 9' and 6 satisfy (2.19), and (3.2) can be written as 

\-Att I A{e',a)v{a)da-q{\)\ = \ j e^-i^' x)<l{x)p{v)dx\ <\\q\\a\\p{v)\\a (3.3) 



where q{X) is defined (2.20) and || • ||q = || • ||i,2(-g^). 

Let us formulate the result: 
Lemma 10: If q G Qa, v & L'^{S'^) and (2.19) holds then (3.3) holds. 

Proof of Theorem A: If one uses (3.3), chooses = iy{a,9) such that (2.17) holds, and uses (2.23), one 
obtains formula (2.21), and (2.22) follows from (2.21). Theorem A is proved. 

2. In the proof of Theorem B the following auxiliary results are important. 
Lemma 11: One has 

\q6 - q{X)\ < MMO',a) - A{0',a)\\^,^^,^a{us) + \\p{iy6)\\a\\q\\a, (3.4) 

\\As{9', a) - A{e', Q!)||z,2(52) < cexp{Kai)iJ,{5) < cexp(K6)/x(^), (3.5) 
where As{e',a) is given by (2.25), A(e',a) is given by (2.10), e'{6) and 9(6) satisfy (2.19), 

K := \Im 0{6)\, i^{S) := ln[--fN{5)], 7 := In ^ > 0, 

aius) := 111^511^2(52), 

q{X) and qs are defined in (2.20) and (2.35), and N{6) is defined by (2.24). 
We collect proofs at the end of this section. 

3. In the proof of Theorems B and C the following results are used. They are of independent interest. 
Lemma 12: If q G Qa, 9 G M, \9\ 1 then there exists tp which satisfies the equation 

iqil) = QinB?, V = exp(i6'-x')[l + i?(x,6')], (3.6) 

and 

\\R\\l-{d) < c\9\-\ \9\ ^ 00 (3.7) 
where D (Z B? is an arbitrary bounded region, the constant c in (3.7) depends on diam D and on \\q\\a '■= 

The following estimate holds [6] : 

ll^ll-(^)<^^- (3-70 
rr 

Let i := iq, Noii) := {w: ew = OmD,w G H'^{D)}, where H'^iD) is the usual Sobolev space, D c 
is an arbitrary bounded region with a smooth boundary. 

Lemma 13: Let w E Nd{£) and e > be an arbitrary small number. Then there exists i'£{a) such that 



\\w - 



j u{x,a)ys{a)da\\L2i^jj) <e (3.8) 



Here u(x, a) is the scattering solution defined in (1.1 ), (1.2). 

Lemma 14: Suppose w = tp and (3.8) holds with w = ip, Im 9 ^ 0. Then 



liml|z.,(a)||=oo, (3.9) 



We\\> cexp{Kd/2), K := \Im 0\, d = diam D (3.10) 



Suppose that 



\\'il}{x,9) — I u{x,a)ve{oi) da\\ <e 



(3.11) 



where the norm in (3.11) is defined in (2.14). This norm is not to be confused with the norm in L'^{S'^). It 
is clear from the context what norm is meant. 
Consider the problem of finding 

iid\\u,{a)\\:=j{e,e) (3.12) 

where the infimum is taken over the functions i^e(a) € L'^{S'^) which satisfy (3.11). 
If 0, € M, is arbitrary large fixed, e > 0, and 



n{e) = exp{(lnlne-i)[l + o(l)]}, e ^ (3.13) 



then the following estimate holds. 
Lemma 15: One has 



j{e, 9) < cexp(Kr) (^ ^"^^^^^ l ^ ^ n''{e), r > b, k = | Jm 0\ (3.14) 



er 



where r > b is arbitrary and c does not depend on e and 6. 
Corollary: Minimizing in r > b one obtains 



j{s, 6) < c(2«)"(^)n2(e), k = \Im 0\ (3.14') 



If N{k) is the asymptotic solution to the equation 

e5«;\^ exp(— «;6) 



, = C2 , K ^ +00 (3.15) 

jy J K 

then, with £(k) := k""'^ exp(— k;6), the following estimate holds. 
Lemma 16: One has 

j{s{K), 6) < ciV3(K)(2K)2^('^)+i, K := I Jm ^| ^ oo, (3.16) 
and N{k) is of the order k as «; — > oo; 

ebK < N{k) < e%K, oo (3.16') 

Thus 

j{e{K), 6) < ck^(2k)2«'^« < ck^(2k)20«'«. (3.16") 

4. Let q G Be, where Be is defined in (1.5). The scattering solution satisfies the equation 

{I + Tg)u = exp{ia- x), TgU := [ — y^q{y)u{y)dy (3.17) 

J 4-K\x-y\ 

Consider the operator I + Tq as an operator in L^{Ba). Note that Tgi L^{Ba) L^{Ba) is compact. 
Lemma 17: The operator I + Tq is an isomorphism of L^{Ba) onto L^{Ba) and 

sup ||(/ + Tg)-i|| <c (3.18) 

qeBc 

where the constant c depends on C and does not depend on q G Be ■ 

5. The following estimate holds (cf. (2.9)): 

er 



\Mr)\<cr-'/'{^^^j j^i^^, 0<r<ro, e>0 (3.19) 



where je{r) is the spherical Bessel function. 

6. Let us prove the lemmas of this section. 
Proof of Lemma 11. One has, using (2.35), 



qs = -4tv j[As{9'{S),a)-A{9'{6),a)]iy5{a)da-4TT j A{e' {5),a)v5{a) da 

By (3.3) one obtains from (3.20) the estimate 

\qs - q{X)\ < Att ( \As{e'{6),a)-A{e'{6),a)\\vs{a)\da+ 



+ Ma\\p{u5)\\a < 4-K\\M0'{5),a) - A{e'{5),a)\\ \\us{a)\\ + |k|U||p(i^^)|U 



This is the estimate (3.4). Let us prove (3.5) 
One has 



N{S) 



\\As{e'{5),a)-A{e'{5),a)\\ < \\ ^ \Asi{a) - A,{a)\ \Ye{e')\ ||+ 



+11 \Mo^)\\ym\\\--=h+h 

l=N(S) + l 

Using estimate (2.8) with r = a\ > a, and (2.9), and taking into account that |Im 9'\ = |Im 
6' - 6 = X, Xe R^, one obtains 



h < cSexp{Kai)N^{d) 



(. 



eai 



\2N{6) + 1 



v/]vp) 



eai 



2W(i5) + l 



< c(5exp(Kai)A^^((5) 



eai/2) 



Here we took into account that there are 



NiS) 

Y{2i+l) = [NiS) + lf 



spherical harmonics with < £ < N{S). 

Furthermore, using (2.7), (2.8) and (2.9) one gets 



E 

e=N{s)+i 



2£ + l 



^ '"""iT^ < cN^S) )"'^xp(.ar) 
^^ + ^f^\— ^ Vary 
V 2^+1 ; 



Combining (3.23) and (3.24) one has, with N = N{6), a/ai := s, < s < 1, eai/2 := t, 



h+h< cexp(«;ai) {6N^[ — ] + N^l'^s 



N 



N 



Consider, for a fixed small 5 > 0, the minimization problem 



5N^ 



N 



+ iVV^s^ = inf 



where the infimum is taken over N, N ^ 1, that is for sufficiently large A''. Write (3.26) as 



Sexp{N\nN-N\nt + 3lnN)+exp{N\ns + -lnN) = inf 



(3.27) 



To find the infimum in N, let us differentiate (3.27) with respect to N and equate the resulting expression 
to zero. This yields 



(5exp(7V lniV-iVlnt + 31n7V) IniV 



l + O 



1 



In TV 



exp(7Vlns+ ilniV)[lns + ^] = (3.28) 



Thus 



^"^ =exp(jVlniV- jVlnt + 31njV- jVlns- ^InTV) ^"^ [^l + ^dnAf 



In - = NlnN 





1 + 



1 



Inln- =ln7V + lnlnAr + 





1 



In TV 



InN 

= (In TV) 



N ^ oo 



1 + 



In In TV 
In TV 



Divide (3.30) by (3.31) to obtain 



N = N{S):=^^[l + o{l)], 5^0 



(3.29) 
(3.30) 
(3.31) 

(3.32) 



Denote by jj,{S) the infimum in (3.27) which is the value of the left-hand side of (3.27) at TV = N{6). Then, 
using (3.27) and (3.28), one gets, with 7 := ln(s^^). 



A((5)=exp{-7TV((5)[l + o(l)]}<^ 1 + 



where n{S) := exp[— 7TV(5)]. 

Prom (3.32), (3.25) and (3.33) one obtains 



2N{5) 



[lnN{6)] [l + o(^) 



< cii{5), 



(3.33) 



/i + /2 < cexp(Avai)/u((5) < cexp(K6)/Ui(5) 



(3.34) 



where 71 := ln[ai /(a\/2)] < 7, 6 > a\. Therefore estimate (3.5) is proved and the proof of Lemma 11 is 
complete. □ 
Proof of Lemma 12. Substitute V of the form (3.6) into the equation (3.6) to get the equation for R: 



Define the operator 



where 



Then 



LR := (V^ + 2ie ■ V)R = q{x)R + q{x) in R^ 
1 f cxp{iX ■ x)f{X) 



(27r) 



A2 + 2A • 61 



dX:= w 



R3 



/(A) := J exp(iA • x)f{x) dx 

Lw = fm R^ 
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(3.35) 



(3.36) 



(3.37) 



(3.38) 



In Lemma 12 the existence of ijj satisfying (3.6) and (3.7) is claimed. Uniqueness of such a ijj is not discussed. 
Therefore, Lemma 12 will be proved if one demonstrates the existence of a solution to the equation 

R = L-^qR + L-^q in (3.39) 

such that R satisfies (3.7). 

Suppose that the following estimate is established 

\\L-'f\\L-^D,)<c\e\-'\\f\\L2(D), eeM, 1^1 »1 (3.40) 

for any / e L'^{R^) vanishing outside of a bounded domain D , D C Ba C Di , Di is a.n arbitrary bounded 
domain in R^, and c depends on D and Di but not on 6 (we will prove estimate (3.40) later: see estimate 
(3.47) below). Then equations (3.39) and (3.40) imply, with || • \\a = \\ ■ h^Ba), 

\\R\U < c\0\-'\\qR\U + c\e\-'\\q\U < c\e\-'\\qhoo^B^) \\R\U + c\e\-'\\q\\a (3.41) 
Therefore, for \9\ ^ 1 such that c|6'|~^||g||L=o(5^) < 1, one has 

\\R\\a<c,\en\q\\a (3.42) 

This implies (3.7) under the additional assumption q G L°°{Ba). 

To complete the proof of Lemma 12 under this additional assumption it is sufficient to prove (3.40). 
This estimate is an immediate consequence of the following result, which can be found in [0, vol. II, pp. 17, 
34]: Let P{d) be a partial differential expression with constant coefficients, d :— —i grad, D C i?", n> 2, D 
is a bounded domain, Go{x) is a regular fundamental solution: P{d)Go = S{x — y), ri{x) G C^{R^), 7] = I 
in a neighborhood of the set {x — y}, x,y G D\, and G\{x — y) = ri{x — y)Go{x — y). Then 



/ 



Gi{x-y)f{y)dy = J Go{x-y)f{y)dy, D c D^, x G Di (3.43) 



D D 



and 



where 



sup P(A)|Gi(A)| < 00 (3.44) 



V{X):={J2 \P^'H^)\Y^' (3.45) 
|j|>o 

and j is a multiindex. If P{d) = i := + 2ie ■ V, then L{X) = + 2A • 61 and 

£(A) = (^ |(A2+2A-0)«nV2 = 

bl>o (3.46) 
={\X^ + 2A • Of + 4|A + Of + 12}V2 > 2|Im e\ 



Therefore 



l|i-Vlli.(z,,) = / 1 / G,{x-y)f{y)dy\^dx<^j d\\G,{\)\Xf{\)? < 

Di D R3 



(3.47) 



Estimate (3.47) is identical with (3.40). Therefore Lemma 12 is proved under the additional assumption 
q G L°°{Ba)- Without this assumption estimate (3.7) is proved in [6] and in [1]. The argument in [6] is more 
complicated. It uses estimate (3.7'): 

\\L-^f\\Loo(D,)<c\e\-^'\\n\e\f'^\\f\\L.^n), 1^1 »1, OgM (3.48) 
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In place of the estimate (3.41) one uses the estimate 

||i?IU~(D.) < c{\9\-^ In 1^1)1/^ lkl|L^(B.)||i?|U~(D,) + c{\e\-^ In \e\f'^ (3.49) 

which impUes 

\\R\\L^(D,)<c{\e\-Hn\e\f'\ qGQ,, 9GM, 1^1 »1 (3.50) 
Using (3.50) and the first inequahty (3.41) one obtains 

\\R\\a < c\e\-'\\R\\Loo^D,)\\q\\a + c\e\-^q\\ < c^\e\-' (3.51) 

Lemma 12 is proved. □ 
Proof of Lemma 13. Suppose that w ^0 and (3.8) is false. Then w is orthogonal in L'^{S^) to all functions 

of the form 

J u{x,a)i/{a) da, 

that is 



D 

Thus 



= J w{x) J u{x,a)iy{a)dadx fov all u G {S'^) (3.52) 
^ ~ J w{x)u{x,a) dx for all a G (3.53) 



It is proved in [3, p. 46] that 

G{x,y) = ^-^p^u{x,a)+o{\y\-^)as\y\^c^, ^ = -a (3.54) 
4TT\y\ \y\ 

where G(x, y) is the resolvent kernel for iq. £gG = —6{x — y) in , G satisfies the radiation condition, Iq is 
defined in (1.1). Therefore (3.53) implies 

= y w{x)G{x, y) dx := h{y) for all y e D' (3.55) 

D 

The function h{y) satisfies the equations 

iqh = -w{x) in R^ (3.56) 

Since h G Hf^^, it follows from (3.55) that 

h = hN = ()ondD (3.57) 

Multiply (3.56) by w{x), integrate over D and then by parts, using (3.57) and the equation £q'w = in D, 
to get 

= y \w\^dx (3.58) 

D 

Therefore ui = 0. This contradiction proves Lemma 13. □ 
Proof of Lemma 14. Suppose (3.9) is false. Then 

||z/e(a)|| < c for all e e (0,eo), £o > (3.59) 
12 



Choose a weakly eonvergent in L'^{S'^) subsequence and denote it v^(a) again, £ — s- 0. Pass to the hmit 
e — » in (3.8) with w = ip{x,9) and this subsequence I'eict), i^eict) '^(a) weakly in LP'{S'^). The result is 

V'(a;, 6)= j u{x, a)v{a) da inD, Im 6i 7^ (3.60) 

Since both sides in (3.60) solve the elliptic equation 

tqip = in i?^ (*) 

they are identical in B? by the unique continuation property for the solution to (*). This is a contradiction 
since 'ip{x,9) grows exponentially as |a;| — > 00 in some directions, while the integral 



u{x,a)iy{a) da\ < c||z^|| 

is bounded. Here we used the well known estimate 

sup \u{x,a)\<c (3.61) 

Estimate (3.9) is proved. In order to prove (3.10) one deduces from (3.8) with w = tp{x,9) that 

II j u{x,a)i's{a,6) da\\L2(D) > ||V'I|l2(£)) - e> cexp(Kd/2), k = jlm e\ (3.62) 

where d =diam D. Assume that for some £ > the estimate (3.10) is false. Then there is a sequence On S M, 
|^„| 00, such that 

||z.e||exp(-K„d/2)^0 (3.63) 
This contradicts inequality (3.62). Indeed, from (3.62) it follows that 

< c < exp(-K„d/2) II j u{x,a)i/e{a,9n)da)\\L2(D) < ciexp{-Knd/2)\\i/s{a,9n)\\L^(S'^) (3.64) 

where ci > does not depend on e or On, it depends on maxr^^ns^cees^ \u{x,a)\ (see (3.61)) and on 
(meas Z?)^/^. Since (3.64) contradicts (3.63), Lemma 14 is proved. □ 

We will use Lemma 17 in the proof of Lemmas 15 and 16. Therefore let us first prove Lemma 17. 
Proof of Lemma 17: First, let us prove that N{I + Tq) = {0}, where N{A) is the null space of a linear 
operator A. It is easy to see that, for q € Qa, the operator Tq in (3.17) is a linear compact operator on 
L'^{Ba). Therefore (I + Tg)^^ exists and is bounded by the Fredholm alternative if N{I + Tq) — {0}. Suppose 
u! + TqW = in Ba- Define w{x) to be —TgW for any x G R^. Then w + TgW = in R^. Therefore w solves 
the problem: 

^ - iwj 0, \x\ 00 (3.65) 

It is well known [19] that (3.65) implies w = 0. Therefore N{I + Tq} = {0} and the operator (7 + Tg)-^ 
exists and is bounded in LP'[Ba). Let us prove (3.18). Assume that (3.18) is false. Then there is a sequence 
qn G Be such that 

\\{I + Tn)-^\\>n, Tn:=Tq^. (*) 

We prove that this is impossible. Indeed, since Be is a compact set in L'^{Ba), one can select a convergent 
in L'^{Ba) subsequence which we denote again by g„, ||g„ — g||o ^ as n — > 00. One can check that 

\\Tq^-Tq^\\<c\\qi-q2\\a (3.66) 
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Indeed 

from which (3.66) follows with c = c(a). Therefore, ii Qn q in L?(Ba) then ||T„ — T|| — > as n — > oo. One 
has, using bounded invertibility of / + T, 

II (7 + T„)-i II = II (7 + T + T„ - T)-i II = II (7 + T)-i [7 + (T„ - T){I + T)-^]-^\\ < c (3.67) 

Here the inequality holds for all sufficiently large n, for example, for such n that ||T„ — T|| ||(7 + T)~^|| < 1. 

Since (3.67) contradicts (*), Lemma 17 is proved. □ 
Proof of Lemma 15: Assume that 9 S M, \9\ ^ 1, is fixed. First, note that one can find v{a) such that 
1 1^1 1 is arbitrary large but 

II j u{x,a)v{a)da\\i,2(jj) <5, (*) 

where 5 > is an arbitrary small number. Indeed, u = Buq, uq := cxp(ia • x) and B = (7 + T)~^ is a 
bounded linear operator in L^(£>), ||B|| + ||B-i|| < c (see Lemma 17). Since \\h\\ = \\B-^Bh\\ < \\B-^\\ 
\\Bh\\, one concludes from (*) that 

c||y ex.p{ia ■ x)i'{a) da\\L2^jj-^ < S, (**) 

where c does not depend on S. Conversely, (**) implies (*) (with ciS in place of 5) since _B is a bounded 
operator. Take i^e{a) = ceYe{a), ||z^|| = Take |q| ^ oo as £ ^ oo so that \ciji{d)\ ^ as i ^ oo where 
d =diam D. Such a choice of q is possible since \je{b)\ — > as £ ^ oo (see (3.19)). So, ||i/£(a)|| — > oo, and 

II j u{x,a)Pi{a)da\\L2(^jj^ <c\\ j exp{ia ■ x)ceYe{a) da\\L2(^D) = c'\ce\\\je{r)\\L2(^D) ^ as £ ^ oo 

The claim is proved. That is why we are looking for ^^(q) which satisfies (3.11) and has minimal norm. 
Let us note that ip = Bip, where (p solves the equation 

\7^ip + (p = OmR^ (3.68) 

Indeed, since (p = B~'^t/j = (7 + T)ip, one has (V^ + l)ip = (V^ + 1)(7 + T)ip = qip - qip = as claimed. 
Inequality (3.11) implies, as in the passage from (*) to (**), that 

||<^ — J exp(iQ: • x)i's{a) da\\ < ce := ei, 6 > ai > aV2, (3.69) 

where c does not depend on s and 6, the norm || • || in this argument is equivalent to the norm || • ||i,2(Bj^), 
and 

^ = V> + TV> = ^+ / ^^^f^^^g(t/)V(2/, 0) dy (3.70) 
J ^■K\x-y\ 

Ba 

Since h > a and \9\ ;» 1, \\ip\\b is of order of which, in turn, is of order of || exp(i0 • a;)||b as 16*1 oo. 
Also, ||<p||&, ll^llb and || cxp(i^ • a;)||6 are of order || exp(i^ • a;)|| as |^| — > oo, ^ e M, where || • || is defined in 
(2.14). Since ip solves (3.68) one can write 

•DC' 

'^(^) =^(^de{r)Ye{a'), r = \x\, a' = x/r (3.71) 
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Also 

oo 

e^p{ia- x) ='^4'Ki^je{r)Yi{a')Yc{a), r = \x\, a' = x/r (3.72) 
Thus (3.69) can be written as 

oo ^ 

\ai - Ani'^iyeel^bl < e\, h\ := / r^\u{r)'^ dr Vei ~ {iye,Ye)L^s^) (3.73) 

One can write an equation similar to (3.17) for tjj: 

iI + T)i; = exp{ie-x), := r,i^ := j G{x,y,0)q^ dy (3.74) 

Here G solves the equation 

{V'^ + 1)G = -S{x), G = exp{i0-x)g (3.75) 
72 , o.-^ V7^„_ „_ ^ f cxp{iX ■ x) 



Lg := (V^ + 210 • V), = -S^x), g=^J ^^^^ dX (3.76) 



The operator (J + r)~^ exists and is bounded in L^{Ba) if \9\ ^ 1. To prove this, it is sufficient, as in the 
proof of Lemma 17, to prove that N{I + F) = {0}. Suppose w + Tw = in Ba- Then w solves equation 
(3.65) and satisfies the following condition at infinity: w = exp{i9 ■ x)v, where 

V = L-^qv in (3.77) 

and is defined in (3.36). The relation (3.77) plays the role of the condition at infinity. Using estimate 
(3.40) one derives from (3.77) restricted to Ba that v = provided that \9\ ^ 1, namely, \9\ is so large that 
\\L~^q\\i^2(^Ba)^L2(Ba) < 1- Thus V = 0, w = 0, and the operator (/ + T)~^ is bounded in LP'{Ba). Since T is 
a compact linear integral operator in L'^{Ba) and / + F is injective, one can write (7 + F)"^ = / + Fi where 
Fi is a compact linear integral operator in L'^^Bg). Therefore 

ij) = exp{i0 ■ x) +Ti exp{i9 ■ x) = exp{i9 ■ x) + J Ti{x, y) exp{i9 ■ y) dy (3.78) 

and from (3.70) one obtains 

ip = exp{i9 • a;) + Fi exp{ie ■ x) +Texp{ie ■ x) + TFi exp(i^ • x) (3.79) 

It follows from (3.79) that exp(z^-a;) is the main term of <^ as \9\ — > oo in the region \x\ > a. Since ip = ip+Ttp 
and IV'I < cexp(K;r) one has \ip\ < cexp{Kr), r > a. Thus 



and 



WfWhis-) = Yl < cexp(2«r), r>a (3.80) 

\ai\ < cexp(/tr)|j£(r)|~\ r > a (3.81) 



Define 

ly^e ■■= {ini^y'^ae for / < N, v^e := for ^ > TV (3.82) 
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Then, using (3.80) and (2.9), one gets 



N 



N 



II ii2 V^i |2 Y^i |2 , exp(2Kr) 2£ + 1 
11^^.11 =}_^Wee\ < 1^ < 



47r ^ |j^(r)|2 



<cexp(2Kr) ^ ( ) (2^ + 1)^ < cexp(2Kr) ( ] 



er 



er 



Let us fix r > 6 and estimate be in (3.73) using (3.19): 

^ 21+1 



(3.83) 



(3.84) 



From (3.81), (3.82) and (3.84) one concludes that (3.73) holds if 



c exp(2K;r) 



t=N+l 



X 2^+1 1 

er \ 1 



2^ + 1 7 2^+1 



eh 1 



< cexp(2Kr) ( 



- 1 < cexp(2K;r) I - ) < £i, r > o 



2^+1 



2N 



Write (3.85) as 



(b\ 
-j < £2, £2 = £i/a/c := C2£, A'^ > 1 

Since r > 6 is arbitrary, one can minimize the left-hand side of (3.86) in r and get 



inf 

r>b 



exp(«;r) 



N' 



= exp{N) 



N 



the infimum being attained at r = N/k. Consider the equality in (3.86) with this r: 



ebK 



N 



= £2, £2^0 



(3.85) 



(3.86) 



(3.87) 



Let us solve (3.87) for N asymptotically as £2 0, for a fixed n. Write (3.87) as 

'In In A" 



lnln£2-i = InTV + InlnTV - + o f ^ ) = (IniV) 



IniV 



IniV 



1 + 



InN 



Therefore 
Thus 



InTV = (lnln£2"i)[l + o(l)], £2^0 
A/' = exp{(lnln£^i)[l + o(l)]}, £2 = C2£, £^0 



(3.88) 



(3.89) 



(3.90) 



With N = N{e) given by (3.90) formula (3.83) yields an upper estimate for and therefore for j{£, 9) for 
arbitrary large fixed 9 G M and £ ^ 0: 



i(£, ^) < cexp(K;r) 



27V + 1 



er 



N 



(3.91) 



16 



where TV is given by (3.90). Minimize the right-hand size of (3.91) in r, r > b. As above, the minimum is 
attained at r = N/k, and (3.91) with this r yields 

/ON + ] \ ^ 

j{e, 0)<c( J K < c(2K)^iV2 (3.92) 



V N 

Since N in (3.90) can be written as n{e) in (3.13), Lemma 15 is proved. □ 
Proof of Lemma 16: If £ = exp{—Kb)/K, k = |Im 9\, then (3.87) takes the form 

e_bny ^ ^^eM-^b) ^^^^^ 

Let N — N{k) solve (3.93) asymptotically as k +x. Denote iy^(a,9) corresponding to e = £(«) := 
K~^exp(— k;6) by i/^. The Fourier coefficients of are given by (3.82) with N = N{k). From (3.83) one 
derives 

Minimizing in r > 6 the right-hand side of (3.94) one obtains, at r = {2N{k) + 1)/(2k), 

exp[2jV(« ) + 1] f 2N{k) + 1 ^ 

{ 2N{k) + 1 \ 

=c exp[27V(K) + 1] ( J (k) = cN^ (k) (2k)2^('^)+i 



(3.95) 

2N{k) + 1 



Estimate (3.16) is proved. To prove (3.16'), note that if N{k) < ebn or N{k,) > e^bn then (3.15) cannot hold 
as «; — > oo. Indeed, as k — > oo then (ebK/N)^ > 1 > C2 exp(— «;6)/k for N{k) < ebn and 



febny fly 

(-^j '^("j < exp(— e^few) < C2K ^ exp(— k;6) for A''(k) > e^^K. 
Lemma 16 is proved. 

IV. Proofs of Lemmas 19 
Proof of Lemma 1: One has, using (3.72) and the orthonormality of the spherical harmonics 



\-^em{oi)\ = j q{x)u{x,a) j exp{-i0 ■ x)Yem{fi) dp dx\ 
< c J \q{x)je{r)Yem{a')\dx, r = \x\, a' = x/r 



(4.1) 



Ba 

and c is the constant from (3.61). From (4.1) using (3.19) one obtains 

1/2 



\AUo^)\ < clklU y Au{T)\^dT\ < era (^^j (4.2) 



Lemma 1 is proved. 

Proof of Lemma 2: Using the formula 



j exp{i0-ar)Ye{a)da = 4TTi'^jiir)Ye{9), OeM (4.3) 
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where r > is arbitrary, one applies the Cauchy inequality to the integral in (4.3) and obtains 

Lemma 2 is proved. □ 
Proof of Lemma 3: The function w, defined by the right-hand side of (2.15), solves the Helmholtz equation 

(V^ + l)w; = for r > a (4.5) 

and has the asymptotics 

w = ex.p{ia ■ x) + A{a' ,a)r~^ ex.p{ir) + o{r~^), i — > oo, a' = x/r (4.6) 

Therefore the function v := w — u{x, a), where u{x, a) is the scattering solution, solves equation (4.5) in the 
region r > a and 

V = o{r~^), r ^ oo (4-6') 

By the Rellich-type uniqueness lemma (see e.g. [3, p. 24]) one has v = for r > a. Lemma 3 is proved. □ 
Remetrk: Lemma 3 is well known. 

Proof of Lemma 4: One has, with if := J u{x, a)i'{a) da, 

S2 



•= II exp(— • x) J u{x, a)i'{a) da — 1|| = 



(4.7) 



S2 

= 11 exp(— 16* ■ x)[(p — ex.p{i6 ■ x) — ex.p{i9 ■ x)R{x, 9) + exp(i6' • x)R{x, 0)]\\ < 
\\eM-iO-x)[^-il:]\\ + \\R{x,e)\\ 

By Lemma 13, with w = ij) and il) as in Lemma 12, one can find v = v{a, 9) such that 

^ exp(-«6) ^ (4.8) 

K 

With this choice of v one has 

||exp(-i6»-a;)[(y9-^]|| < exp(K&)||(^ - -i/;!! < k"^ (4.9) 

Since 6* G M one has \9\/k ^ \/2 as ^ ^ 00. Prom this, (4.7), (4.8) and (4.9) the estimate (2.16) follows. 
Lemma 4 is proved. □ 
Proof of Lemma 5: Let us choose the coordinate system so that A = te^, where t> Q and Cj, j = 1, 2, 3, is 
an orthonormal basis of i?^ . Let 9' := ^te^+v, 9 := —^te^+v, where w-w = 1— 1^/4, W3 = 0, \v\ — !■ 00, v S C^. 
Clearly there are infinitely many such v. li t 1, then one can choose 9' and 9 in M' in the above proof 
Lemma 5 is proved. □ 

The reader can also consult [4] for another proof. 
The proof of Lemma 6 is given in Section V. 

Let us prove the remaining lemmas. 
Proof of Lemma 7: One has 



l|P5(^)ll •= II exp(— • x) J us{x, a)i'{a) da — 1|| 



< II exp(— • x) J u{x, a)u{a) rfa — 1|| + || exp(— • x) J [us{x, a) — u{x, a)]i'{a) da\\ < (4-10) 
— Ilp(^)ll + cexp{Kb)a{i') sup ||u5(a;,a) — ^(x, a)|| 

aeS^ 
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Let us prove that 

sup \\us{x,a) — u{x,a)\\ < cfi{6) (4-11) 

where iJ,{S) is defined in (2.28) and N{5) is defined in (2.24). One has, with A'' = N{5), r = \x\ and a' = x/r, 

N oo 

\\us{x,a)-u{x,a)\\ < Wj^i^sda) - Ae{a)]Ye{a')he{r)\\ + || ^ Ae{a)Ye{a')he{r)\\ := h + h (4.12) 
e=o e=N+i 

Using Parseval's equahty and the assumption (1.3) one gets 

II = m&xJ2\^se{a) - Ae{a)\^ j r'^\hi{r)\^ dr < 4^6'^ / r^\he{r)\^ dr. (4.13) 

Furthermore 

J| < max £ \A,{a)\' [ r'Mr)\' dr (4.14) 



e=N+i 



Define 



One can prove (see [5]) that 



b 

2 



He := / r'\he{r)\Ur (4.15) 



He<c('-^Yiy^ a. = ^ (4.16) 
\ ea2 J V2 



In fact, a stronger estimate holds: 

„ , o.,. _i /2£+]^ 



ff^ < ca?[(2^ + l)£]-i . (4.16') 



eai 



This follows from the known asymptotics of Hankel's functions: he{r) ^ —i{rt) "'^ (^^7^)^^°^ as ^ — > 00. 
Prom (4.13)-(4.16') and (2.7) one gets 



h+h<c 



27V + l\^+i /a^^' 



(4.17) 



Minimizing with respect to N the expression in brackets in (4.17) for a small fixed (5 > one obtains that 

the minimum is c/ii((5) with /ii((5) defined in (2.28) and N = N{5) defined in (2.24) (cf. [5]). Inequality 

(2.30) is proved. Inequality (2.31) can be proved similarly. Lemma 7 is proved. □ 
Proof of Lemma 8: Use (2.30) to conclude that (2.32) holds if 

||p(i.)|| + a{v) exp(K6)^(<5) < c\e\-^ (4.18) 

Using Lemma 13 with e = \6\~^ exp(fi;6), k = |Im 9\ and D = B},, choose u = y{a, 6) such that 

\\^ - i'W < \0\~^ exp(-K&) (4.19) 

where 

ip := u{x,a)iy{a,0)da (4.20) 
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Then, using (3.6), (3.7) and (4.19), one gets 



||p(i/)|| =||cxp(-«0-a:)^-l|| < \\cM-iO-x)[ip-ij]\\ + \\R\\ < 

< exp(/t6)||^ - VII + c\e\-^ < \e\-\c + 1) 



(4.21) 



Choose 9 = 6{S) such that 

\e\a{iy{a, 6)) exp(K6) = l//x(5) (4.22) 
Note that l-t~^{S) — > oo as ^ — > 0. We claim that 

a(i/(a,6')) ^ oo as |6i| ^ oo, 6* e M (4.23) 

Prom (4.23) it follows that equation (4.22) has a solution 6{S) such that 

\e{6)\ -^ooasS^O (4.24) 

To finish the proof, let us check that (4.23) holds. Assume the contrary, that is 

a{u{a,0)) <c' (4.25) 

where c' docs not depend on 6. Using (4.25) choose a weakly convergent in L'^{S'^) subsequence which is 
denoted by i^(q, ^„):=f„, 

v{a,en) i^{a) weakly in L^{S^) as \0n\ oo (4.26) 
Pass to the limit in (4.19) to get 

\\7p{x,9n)-ip{x,0n)\\^O, n^oo (4.27) 

where 

ip{x,6n) '■= J u{x,a)i'{a,6n) da (4.28) 
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This is a contradiction since \\(p\\ < c||t'„||(47r)^/^ < cc'(47r)'^/^ < Ci, while ||V(a;,^n)|| ^ 00 as n — > 00. In 
this estimate c is the constant from (3.61). Lemma 8 is proved. □ 
Remark: The rate of growth of a(j^(a, 6*)) as \9\ — !■ 00, 9 £ M, can be estimated from below by formula 
(3.10) and from above, for some choice of v{a,9), by (3.16)-(3.16"). This allows one to get the following 
lower estimate of 0((5) which gives a refinement of (2.33): 

To prove (4.29) let us solve (4.22), asymptotically as 5 ^ 0, for \9{5)\ using the upper estimate (3.16") for 
a{v,9). Equation (4.22) becomes 

1^1 exp(K6)K^(2K)20«'« = (4.30) 

Note that in order to obtain a lower estimate for \9{8)\ one has to use an upper estimate for a{v) in (4.22). 
Since \9\/k \/2 as \9\ 00, 9 € M, one can write (4.30) as 

exp[K6 + 51nK + 206k In k + 206k; In 2] = exp[7A''((5)] (4.31) 

where N{6) is given in (2.24). Note that the main term of the asymptotic solution to (4.31) will be the same 
for the equation (4.31) with cexp['yN{S)] in place of exp['yN{d)], c =const. Prom (4.31) one derives 



206K(lnK) 



1 + 



1 

luK 



: ^N{6) (4.32) 
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Thus 

«lnK = 7iiV(5)[l + o(l)], <5^0, 11 ■■=^ (4.33) 
Taking In of (4.33) and denoting 'yiN{6) := n yields 

InK + lnlnK = lnn + o(l) (4.34) 
Let us look for the asymptotic solution to (4.34) of the form 

Tl 

K=- — (1 + t), t = o(l)asn^oo (4.35) 
In n 



Substitute (4.35) into (4.34) to get 

Inlnn + f + 0(i2) 



lnn-lnlnn + t + 0(i^)+ln<( (Inn) 

Thus 



1 - 



Inn 



= lnn + o(l) 



t + Oit^)-'^^^^^^l±^=oil) (4.36) 
Inn 

Equation (4.36) implies that t = o(l) as n ^ oo. Therefore formula (4.35) gives an asymptotic solution to 
equation (4.33) as 5 ^ 0. Let us write this solution for references in terms of \0{S)\ = \/2k;(^)[1 + o(l)] as 

where N{5) = \ ln5|/ln | ln5|. Therefore 



Let us formulate the result: 
Lemma 8': One has 

V2ji\\nS\ 
{\n\\nS\) 

where 

\n[ai/a 



e{5) > \e{S)\ = + o{l)] as 5 ^ (4.39) 



71 = > (4.40) 

Proof of Lemma 9: The starting point is the standard resolvent identity 

Gi{x,y)-G2{x,y) = - ( G^{x,z)p{z)G2{z,v)dz (4.41) 



where p := qi — and Gj is the resolvent kernel of the operator £j := V'^ + 1 — qj{x). Let \y\ oo, y/\y\ = P 
in (4.41). Using (3.54) one gets 

ui{x, -13) - U2{x, -P) = - J Gi{x, z)p{z)u2{z, -13) dz (4.42) 

Let |a;| oo, x/\x\ = —a in (4.42) and use (1.2) and (3.54) again to get 

A^i-a, -13) - A2{-a, -(3) = ^ j p(z)ui(z, a)u2{z, -(3) dz (4.43) 
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By the well known reciprocity property Aj(—a, —ft) = Aj{P,a), so (4.43) is identical with (2.40). The 
original proof of Lemma 9, given in [9], was a little longer. □ 

V. Proof of Lemma 6 
This proof requires some preparation. Consider the equation 



Let 



Lp := (V^ + 2i6» • V)p = w in i^^ peC^(Sr), 6' S M (5.1 

^(0 := + 2/3 • <e, P-=he h:=\e\-^ (5.3 

(}-(} = h\ \(3\ = 1 (5.4; 

N:={^:P{0=0, ^GR'} (5.5 

Nh := {C: ^ e i^^ dist (C, A^) < h} (5.6 

at;; := \ Nh (5.?: 

Let f3 = m + ijj., m, e i?^. Then 

iV = {^:|^ + m| = |m|, /i • ^ = 0, ^ e i?^} (5.8 
Thus iV is a circle. Let P(0 = ^^1(6 + i-f2(C), where Pi(0 :=Re P(^). Note that 

dPi(^)5^0for^e (5.9 
where dPi is the differential of Pi. Equation (5.1) can be written as 

P{hD)p := [{hDf + 2/3 • hD]p = -h\ (5.10 



Define 



FhU := u{^) := (277/1)-^/^ j ^j-j.) exp(-i^ • xh''^) dx (5.11 

M(a;) = (27rft,)-3/2 j gxp(i^ • x/i'^) (5.12 

ihd^.u{^) = XjU (5.13 

P^-i/ia,«(a;)} = ^ju(0 (5.14 

Let us denote in this section 

IIpII := l|p||L^(fl3) (5.15 

\\p\Wm ■■= \\p\\L-^iB,,\B,,), < 61 < 62 (5.16 

\\g{hD)p\\ := uViTemow (5.17: 
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First, wc need 

Lemma 18: Any solution p G CQ{Br) to (5.10) satisfies the inequality 



h\\{hDfp\\ < c\\P{hD)p\\ V/i e (0, ho) (5.18) 

where ho > is a fixed sufEciently small number. 
Proof: Inequality (5.18) can be written as 

h\\ii+e)p\\<4Pm\\ (5.19) 

where the definition (5.17) and Parseval's equality are used. 
Let ^&N'f^. Then 

h{i+e)<c\pm (5.20) 

so that 

h-" j{i+ef\p?d^<c^ j \pm'\m\'di<c' j \p{om\'d^ = c'j \p{hD)p\^dx (5.21) 

If ^ e Nfi, then one introduces local coordinates in which the set N has the equations 

t = 0, 6=0, i = Pi(e), (5.22) 

the ^-axis being chosen along the vector p, defined by the equation (3 = m + ipi. These local coordinates can 
be introduced because of the condition (5.9). Define 

f--=pmm (5.23) 

Then / = at i = and / e C°°{R^). Let us use the inequality 

h h 

j t-'\m\'dt<Aj \rit)\'dt (5.24) 

-h -h 

proved in Lemma 19 below. This inequality holds for 

/ e C\-h, h), h = const > 0, /(O) = (5.25) 
Applying (5.24) to (5.23) and integrating (5.24) in the remaining variables ^ one gets 

j \m?d^ < c I \ydPiiom]\^d^ < c J lydPiUMiWd^ (5-26) 

Note that the set Nh is bounded. Therefore 

1(1+ ef\m\' < ch' I \m\' d^ (5.27) 

Using Parseval's equality, the assumption supp p C Br, and the S.Bernstein's inequality for derivatives of 
entire functions of exponential type, one gets 

\^dPiiOm\'d^^r' I \P,{i)p\Hx<r^ j \Pm' dx = j \P{hD)p\Ux (5.28) 
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Prom (5.26)-(5.28) one obtains 

/i2 y" (1 + e')'|/5(e)|' rfC < cr^ J \PihD)p\^ dx (5.29) 



Prom (5.21) and (5.29) one gets 



J {l+^^f\p{^)fd^<c^j \P{hD)p\^dx (5.30) 



This inequality implies (5.18). Lemma 18 is proved. 

Lemma 19: Under the assumptions (5.25) inequality (5.24) holds. 

Proof: Inequality (5.24) is similar to the well known Hardy's inequality in which the integration is taken 
over (0, oo). To prove (5.24) one starts with an obvious inequality in which A is an arbitrary real number: 



h 

If - Xt-^f{t)fdt:= AX^ - XB + C (5.31) 

-h 



where 



and 



One has 



h 



A:= Jt-^\ffdt, B:= Jt-'j^\fit)fdt (5.32) 



h 

C := I \ffdt (5.33) 

-h 



B = t-'\f{t)f\>lh + J t-^\f?dt= l-^^^^' +A>A (5.34) 

-h 

Prom (5.31) it follows that 

B"^ < AAC (5.35) 

By (5.34) one Yias B > A and (5.35) implies 

A<AC (5.36) 

This is inequality (5.24). Lemma 19 is proved. □ 

We need one more lemma: 
Lemma 20: Let 

P{hD)p = in Ai (5.37) 

where Ai is a bounded region with a smooth boundary. Let A C Ai, rj{x) S Cq°(Ai), < rf{x) < 1, ri{x) = 1 
in A, wiiere A is a strictly inner subdomain of Ai. Then 

h\\DpU < c\\p\\a, (5.38) 

where \\p\\a := 

Proof: Multiply (5.37) by rjp, take the real part and integrate by parts to get 



h J r]\S/p\'^ dx = J{pVp + pVp)V'ndx + 2Re Uf^j J PjPVdx 
= y IpI^ V^ry dx + 2Re I iPj J pjprj dx 



(5.39) 
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where the summation over the repeated indices is understood. Using the inequalities 

|2ftP|<^|p,f + 2/.-i|p|2 



(5.40) 



\^^v\<c, \Pj\<l (5.41) 
one obtains from (5.39) the following inequality 

h j p\'^ dx < ch j \p\'^dx+^ j r]\pj\'^ dx + 2h-'^ J v\pf dx (5.42) 

Ai Ai Ai Ai 

It follows from (5.42) that 

j \Vpf < ci/i^ j dx + C2 j 'n\p? dx<Ci j \p\^ dx (5.43) 

A Ai Ai Ai 

This inequality implies (5.38). Lemma 20 is proved. □ 
We are now ready to prove Lemma 6. 

Proof of Lemma 6: The equation for the function p defined in (2.12) is 

Lp = qp + qm (5.44) 

where L is defined in (5.1) Write (5.44) as 

P{hD)p = -h'^{qp + q), h := \e\-'^ (5.45) 
Let T] e C^{Bb), < r]{x) < 1, 7]{x) = 1 in Ba,. Clearly 

P{r]p) = {Pr] - 'nP)p - h^r]{qp + q) (5.46) 

Applying (5.18) to (5.46) yields 

h\\{hDf{pr,)\\ < c\\{Pr,-riP)p\\+ch^qh^^sJ\p\\a + ch'\\q\\a (5.47) 
Since r]{x) = 1 in Ba, one obtains 

h\\p\l < h\\{hDf{r]p)\\ < ch^p\l + ch^ + c\\iPr]-vP)p\\ (5.48) 

Thus 

\\p\\a<ch + ch-'\\{Pri-rjP)p\\ (5.49) 

Using the equation 

Dr] = Om Bai (5.50) 

one obtains 

\\{PTj-rjP)p\\ = \\p{hDfri + 2h^Dri ■ Dp + 2h(3p ■ Dt^W < c{h^ + h)\\p\\a,.b + ch^Dp\\a„t (5.51) 

The function p solves equation (5.44) and q{x) = in i?b \ Ba^. Therefore Lemma 20 is applicable and the 
estimate (5.38) yields 

h\\Dp\l,^, < c\\pU-,,b+e (5.52) 
where £ > is an arbitrary small number. 
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Prom (5.49), (5.51)and (5.52) one obtains 

\\p\\a <ch + c||p||ai-e,6+£ (5.53) 

Since £ > is arbitrarily small, one chooses e so small that a < ai — e, where ai > a, and (5.53) implies 
(2.23). Lemma 6 is proved. □ 
The arguments in this section are close to those in [8] . 

VI. Proof of Theorems B and C 
We have already proved Theorem A in Section III.l. 
Proof of Theorem B: Prom (2.35) with 6 = 6(6) and 6' = 9' (5), where 9(6), 6' (6) and us{a) are defined in 
Lemma 8, one has 

qs = -47r y [As{9', a) - A{e' , a)]iysia) da-A-n j A{e' , a)vs{a) da := h + h (6.1) 

The last term in (6.1) is transformed as in (3.2) and (3.3) and one obtains the estimate 

\-At, j A{0',a)vs{a)da-q{X)\ < \\q\\a\\p{'^s)\\a < c\9\-^ (6.2) 

where formula (2.32) was used. 

The first term is estimated by (3.5) and (2.32): 

|Ji| < ca{us)exp{Kb)i^{6) < c\9\-'^ (6.3) 

Prom (6.2) and (6.3) the estimate (2.36) follows. 

Let us prove (2.37). To this end one has to estimate the function &{6), defined in (2.32), from below. 
Such an estimate is given in (4.39) and yields (2.37). Finally, Lemma 17 shows that the constants c in our 
estimates can be chosen uniformly for q € Be- Theorem B is proved. □ 
Proof of Theorem C: The starting point is formula (2.40). Multiply (2.40) by iyi{a,e) and U2{-a',9') 
and integrate over S'^ x S'^ to get 

-4:TT J j A{a' ,a)ux{a,9)v2{-a' ,9')dada' = j p{x)(pi{x,9)<f2{x,9') dx (6.4) 

Here 9, 9' € M, \9\ ^ 1, 9' + 9 = X, X € R^, 

ipj{x,9j) := J Uj{x,a)i'j{a,9j) da, j = 1,2, 9i:=9, 92 ~ 9' (6.5) 

and fj = Vj{a, 9j) are chosen so that 

\\p{'^j)\\ < c\9\-\ p{vj) = exp(-*0,- • x)^, - 1 (6.6) 

where \9i\ |^2|~^ — > 1, as — > oo, and the norm is defined in (2.14). Then, using (6.6) and (2.23), one 
obtains 

j pix)(piip2 dx^ j p{x) exp(iA • x)(l + pi){l + pa) dx = p{X) + 0(16*1)-^) (6.7) 

S„ Ba 

Note that for |A| < Aq one can choose a constant c independent of A, such that 

0{\9\-') < c\9\-\ 1^ ^ 1 as l^il ^ oo (6.8) 
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The left-hand side of (6.4) can be estimated by the Cauchy inequality. Using (2.41) one gets 

|-47ry" J A{a,a)piiy2dada'\<16TT'^Sa{ui)a{u2) (6.9) 

From (6.4) and (6.7)-(6.9) one obtains 

|p(A)| <c[5a(i/i)a(i/2) + |er'] (6.10) 

Let us estimate a{vj). If one takes e = «;~^ exp(— «;6) in (3.11), then estimates (6.6) hold (see formulas 
(4.7)-(4.9)). Therefore, by (3.16"), 

a{uj) < ck^(2k)20*»^, k = |Im = |Im (6.11) 

Thus (6.10) and (6.11) yield 

|p(A) I < c[5k^ (2k)4°'"' + K- 1] (6.12) 
For a small fixed 5 > find the infimum of the right-hand side of (6.12) 

inf \5k^{2k)^^'' + K-^\ := r](6) (6.13) 

k;>0 

One has 

dh{K) + K~^ := (5exp(406KlnK + 406Kln2 + 81nK) = inf (6.14) 
Taking the derivative with respect to k yields, at the point of minimum, 

Sh{K)[minK + 0{l)]= K-"^, K>1 (6.15) 

Let us solve (6.15) asymptotically, as (5 — > 0, for k = k{6). This allows us to estimate ry(^) as 5 — > 0. Taking 
In of (6.15) yields 

In^ =2\nK + mK\nK + mKln2 + 8lnK + ln[minK + 0{l)] = 406K(lnK)[l + o(l)], k^O (6.16) 
d 

This is an equation similar to (4.33). Using the argument given for the proof of Lemma 8', one obtains the 

asymptotic solution (6.16). The role of N{S) is played by ln{l/S), and the role of 72 is played by 1/(406). 
Thus, (4.35) yields 



Equation (6.17) can be written as 



From (6.13)-(6.15) and (6.18) one obtains 

ri{S) = -{1 + o(l)] = 406i^^j^[l + o(l)] as (5 ^ (6.19) 

K I In ^ I 

This and (6.12) yield (2.42). The constant c by Lemma 17 can be chosen uniformly for q G Be- Theorem C 
is proved. □ 
Theorem C is a refinement and an improvement over the result in [14] (see also [31]). 

VII. Summary of the Numerical Procedure, Additional Results and Remarks 
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In this last section wc disciiss the following items: a) stability of the rec;overy of q{x) given q{X) with 
some error, b) discussion of the steps in a possible numerical implementation of the methods for solving the 
ISP developed in this paper, and c) an alternative to (2.32) optimization problem which leads algorithmically 
to a stable approximation of q{X) from the knowledge of noisy data. 

1. Suppose 95 (A) is known such that 

sup MX) - q{X)\ < 6 (7.1) 

|A|<Ao 

where Aq > is a large number, q{X) is defined in (2.20), and a priori it is assumed that q{x) € Q and 

|9(A)| <ci(1 + A2)-^ J>|, XeR" (7.2) 
Given ^^(A) for |A| < Ao, and the numbers ci, d and 6, one wants to estimate q{x) stably as 5 ^ 0. Define 

qs{x) := {2^)^'^ J qs{X)exp{iX- x)dX (7.3) 

|A|<Ao 



Then 



95(2;) = ■^^"p' j i^sW - qW]exp{iX- x)dX- J q{X) exp{iX ■ x) dX+ 

|A|<Ao fl^\BAo 



(2^ / ^'^^^ ^^Pi^^ • ^) '■= h+l2 + q{x) 



(7.4) 



One has 



|/i|<<5A^^ (7.5) 

CXD CSO 

4^ci f X^dX ci [.-2d+2.._ A^ 



Prom (7.4)-(7.6) one gets 



\qs{x) - q{x)\ < ^ + -~ ^ (7.7) 



For a fixed small S one can minimize the right-hand side of (7.7) in Aq and find Xo{S) at which the minimum 
is attained: 



Xo{S)={^y' (7.8) 



This equation gives a practical estimate of the region on which 55 (A) should be known for a stable recovery 
of q{x). With Ao = Ao(^) given in (7.8) one obtains from (7.7) the following estimate 

\qs{x)-q{x)\<coS^-^ (7.9) 

where 



Co := 



1 

+ 



67r2 27r2(2j-3) 



(7.10) 



Let us formulate the result: 

Lemma 21: Let (7.1) and (7.2) hold, the numbers S, C\ and d he known, Xq{5) he defined in (7.8), and 
qs{x) he given in (7.3). Then the error estimate of the recovery of q{x) by the formula (7.3) is given hy (7.9) 
and (7.10). 
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2. Let us discuss the steps in a numerical implementation of the methods for solving the ISPg developed 
in this paper. 

Step 1: Given As{a',a) one calculates As{6\a) by formula (2.25), then us and P5{y) by formulas (2.26), 
(2.27) with e' and 9 satisfying (2.19). 

Step 2: One solves the variational problem (2.32) taking 9 of order given in (2.37). This can be done by 
minimizing the functional \\ps(i')\\. For example, one can look for v of the form 



u = Y,veYe(a) (7.11) 



1=0 

and find the coefEcients from the linear system which one gets from the condition 

\\ps{i^)\\ =mm (7.12) 

If us{a, 6) is an approximate solution to (7.12), one checks if the inequality (2.32) holds with some constant 
c independent of 5. This is done by solving several problems with |, |, f in place of 5. If the inequality 
(2.32) does not hold, one should increase the number n in (7.11) and try to decrease \9\. 

One may do numerical experiments with some q(x), for example q(x) = 1 in Ba, q(x) = outside Ba, 
and get an idea about the values of 9{5) and vs(a) needed for recovery. If formulas (2.38)-(2.39') are used 
then only discrete noisy data As(a'p, oij) are used for recovery. 

3. Let us formulate in conclusion an alternative to (2.32) optimization method for solving the ISP with 
noisy data. The alternative method consists of the following steps: 

1) First, solve the problem 

\\ps(v)\\=mi, ueL\S^) (7.13) 
Denote by d(6,6) the infimum in (7.13) and find ve^sia.) such that 

\\ps{'^e,6)\\<d{S,9) + \9\-' (7.14) 

2) Secondly, solve the problem 

H(6, 6) := inf := ij(S), 6» e M (7.15) 

where 

H(6, 9) := \\ps{i^0,6)\\ + a(iy0,s) exp(Kb)p(S), K = lm \9\, (7.16) 

and find 9i(d) such that 

H{6,9i(6)) <io(6)+S (7.17) 

Define qis by the formula 

gia = -47ry A(6'i((5),a)z/i5(a)da, i^is ■= i^ei{5),s (7.18) 

where O'^ (S) — 6i (S) = A. Then the following error estimate holds 

sup \qu - q{X)\ < c[\ei(S)\-' + 5 + \e(S)\-'] (7.19) 

|A|<Ao 

Here 9(6) satisfies estimate (2.37) and 9i(S) is calculated numerically (see (7.16) and (7.17)). 
Let us prove (7.19). First, note that 

11^(1^15)11 < c[u!(d) +6], c = const > 1 (7.20) 

Indeed, by (2.31), 

Wpii^is)]] < \\psii^i5)\\+caii^i5)exp(Kib)ii(5) < c[u;(5) + S], c > 1, m = |Im 6»i| (7.21) 
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Secondly, as in (6.1)-(6.3), one gets, using Lemma 6, and (3.4), 



1^15 - 9(A)| < c[\\p{vu)\\a + a{vu) exp(Ki&)Ai(<5)] < c[iv{5) + 5 + |0i((5)r^] (7.22) 
Thirdly, let us prove that 

w(^) < c\e{5)\-^ (7.23) 
where 9{5) is defined in (2.32) and satisfies (2.37). One has, using (7.14), 

ij{5) < mi^[d{S,0) + \e\-^ +a{iy0,s)exp{Kb)iJ.{6)] < c\e{6)\-^ (7.24) 

where 0{6) satisfies (2.37). 

From (7.20), (7.22) and (7.23) estimate (7.19) follows with \e{S)\ satisfying (2.37). Let us summarize 
the result assuming that |^i(i5)| ^ oo as ^ ^ 0: 

Lemma 22: The function (7.18) is a stable approximation of q{X) and the error estimate is given in (7.19) 

with \e{S)\ satisfying incquaUty (2.37) and \ei(S)\ defined by (7.17). 
The result of Lemma 22 is an improvement over the result in [18]. 

4. The methods for solving the ISP5 developed in this paper are applicable in many fields. For example, 

the applications to geophysical inverse problems arc discussed in [1,6,20,26], to Maxwell's equations in [1], 
[21], to hyperbolic equations in [1], [22] [30], to inverse conductivity problem in [1], [23-25], to inverse spectral 
problem in [1],[29] and property C , the basis of our theory, was introduced in [27], [6], [28] and [1]. 

5. Finally, let us point out that for a study of ISP s it is necessary to assume that q(x) vanishes outside 
some hall. Indeed, if q(x) G L^{B?) then the function 

r, (^\ - j li^)' if l^^l > a; 

contributes to the scattering amplitude 

A{a' ,a) = —(4:Tr)^^ / exp(—ia' ■ x)u(x,a)q{x) dx 



JJ3 

the quantity which does not exceed 

(47r)~"^c J \q{x)\ da; — > as a — > 00, 

|x|>a 

where c is the constant in (3.61). This contribution becomes indistinguishable from the noise when 

y \q{x)\ dx < 4nc-^6 (7.25) 

|x|>a 

Thus, one cannot recover q(x) in the region |a;| > a, where a is determined by (7.25). 

From (7.25) one can estimate the order of the radius a of the ball in which one can recover the potential 
given the noisy data with the noise level 5. For example, if one knows a priori that \q{x)\ < ci|x|~'' , d' > 3, 
\x\ > 1, then (7.25) implies a^'''' < ((d' - 3)/cic)S, so that a > (d' - 3){cic)-^S^^(-'^'-^\ 
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